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ABSTRACT We solve the Cauchy problem of the Ward model in light-cone coor- 
dinates using the inverse spectral (scattering) method. In particular we show that the 
solution can be constructed by solving a 2 x 2 local matrix Riemann-Hilbert problem 
which is uniquely defined in terms of the initial data. These results are also directly 
applicable to the 2 -f 1 Chiral model. 
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1 INTRODUCTION 

We study the Cauchy problem for the Ward model in light-cone coordinates: 

Qxt Qyy [Qy 

], x,y e R, t >0, (1) 
Q{x,y,0) = Qoix,y), (2) 

where [ , ] denotes the usual matrix commutator, Q{x,y,t) is a traceless 2x2 anti- 
Hermitian matrix and Qo{x,y) is a 2 x 2 anti-Hermitian traceless matrix decaying suffi- 
ciently fast a.s x"^ + y"^ ^ oo. 

We shall solve this problem using the so-called inverse spectral (scattering) method. 
This method is based on the fact that equation (|1|) is the compatibility condition of the 
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following Lax pair, 



/iy - kjj,^ - Q^j^ = 0, (3) 
fit - /cVx - {kQ^ + Qy)fi = 0, keC, (4) 

where ii{x,y,t, /c) is a 2 x 2 matrix. The transformation 

t-y t+y 

x,y,t — > — X, — (5) 

maps equation (|l]) to the Ward model [jl| in laboratory coordinates. The Cauchy problem 
in this model is defined by 

Qtt - Qxx - Qyy = [Qx, Qt] - [Qx, Qy], X,y & R, t> 0, (6) 

Qix,y,0) = Qiix,y), Qt{x,y,0) = Q2ix,y). (7) 

This problem can be solved by using the fact that equation possess the following Lax 
pair 



{k -^)H:,,-2ny- {Q^ + '^^—j^) H = 0, 



2/^t-(^ + ^K-(-gx + ^^y^)/i = 0, keC. (9) 
The Cauchy problem (||), (J^), was studied in p using the Lax pair (H), (H). 

Here we study the Cauchy problem (P, (|^), using the Lax pair (|]), @. We also make 
some remarks about the Cauchy problem (|), 

We note that the transformations 

Qy = —J ^Jf, Qx = —J ^Jy, (10) 

and 

Qx = -{J~'Jt + J~'Jy), Qt - Qy = -J-'Jx, (11) 

map equations (|lD and (|) to equations 

(J-V,)^ = (J-J,)^, (12) 

and 

{j-'J^)x = {j-'Ji^-y)),-yy (13) 
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respectively. Thus, our results are directly applicable to the solutions of the Cauchy prob- 
lem for equations (0) and (|13|). These equations are the 2+1 integrable chiral equations 
in light-cone and laboratory coordinates, respectively. 



In order to simplify the rigorous aspects of our formalism we first assume that Qo{x, y) 
is a Schwartz function which is small in the following sense 

/ Qo{^,y) d^dy-^l, (14) 

where Qo is the fourier transformation of Qq in the x variable. This assumption excludes 
soliton solutions. We then indicate how the formalism can be extended in the case that 
the above assumption is violated. In the case that Qq is sufficient small, the inverse 
spectral method yields a solution of the Ward model in light-cone coordinates through 
the following construction. 



Theorem 1 Let QQ{x,y) , x,y E R be a 2 x 2 anti-Hermitian traceless matrix which is a 
Schwartz function and which is small in the sense of equation ^^}. 

(i) Given Qq{x, y), define /i^(x, y,k), k E = {k E C : Imk > 0} and /i~(x, y, k), 
k G = {k E C : Imk < 0} as the 2x2 matrix valued functions which are the unique 
solutions of the linear integral equations 

/i+(x, y, k)=I+^( rdprdy'-f'dprdy') [ dx' e'^^^^^'+^^y^y'^^QoA^' , y')^i^{x\ y\ k), 
47r\Jo J-oo J-oo Jy JJr 

(15) 

and 

/.-(x, y, k)=I+U fdpf'dy'- rdpHdy') f dx'e'^^'-^'+"~y-y">^QoAx' , y')f^-{x', y', k), 

Att\J-oo J-oo Jo Jy JJr 

(16) 

where I denotes the 2x2 unit matrix. 

(a) Given fi^ define the 2x2 matrix S{x + ky, k), x,y,k E R, by 

I -S = (l -— C dpe'P^''+''y'> I dx' dy' e-*P(^'+'=2'') Q^^^ix' , y') fi+(x', y', k)] x 

V 47r J-oo Ji?2 / 

(l r dpe'P^''+'y^ f dx'dy'e-'P'^^'+'y'^Qo.'{x',y') fM-{x',y',k)) . (17) 

V 47r JO Jij2 / 

(iv) Given S{x + ky,k) define the sectionally holomorphic function M{x,y,t, k) = 
M"*" (x, y, t, k) for k G , M(x, y, t, k) = M~ (x, y, t, k) for k G as the unique solution 
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of the following 2x2 Riemann-Hilbert problem 



M-{x,y,t,k) = M+{x,y,t,k) (^I - S{x + ky + kH,k)^ , k e R, (18) 
detM = l, (19) 



M = I + 0{l/k), k^oo. 



(20) 



(v) Given M{x,y,t, k) define Q as 



Qx{x,y,t) = — [ 



dk M~^{x, y, t, k)S{x, y, t, k). 



(21) 



Then Q solves equation and Q{x,y,0) = Qo{x,y). 

The rest of the paper is organized as follows. In section 2 we derive the main theorem. 
In section 3 we show how the relevant formalism can be extended to include soliton 
solutions. In section 4 we briefly discuss the Cauchy problem for the Ward model in 
laboratory coordinates. 

We now make some remarks about related work. A method for solving the Cauchy 
problem for decaying initial data for integrable evolution equations in one spatial variable 
was discovered in 1967 0. This method which we refer to as the inverse spectral method, 
reduces the solution of the Cauchy problem to the solution of an inverse scattering problem 
for an associated linear eigenvalue equation (namely for the x-part of the associated 
Lax pair). Such an integrable evolution equation in one spatial dimension is the chiral 
equation; the associated x-part of the Lax pair is 



where the eigenfunction fi{x, t, fc) is a 2 x 2 matrix, k is the spectral parameter and Q{x, t) 
is a solution of the chiral equation. 

Each integrable evolution equation in one spatial dimension has several two spatial 
dimensional integrable generalizations. An integrable generalization of the chiral equation 
is (0). A method for solving the Cauchy problem for decaying initial data for integrable 
evolution equations in two spatial variables appeared in the early 1980 (see reviews 0,0). 
For some equations such as the Kadomtsev-Petviashvili I equation, this method is based 




(22) 
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on a nonlocal Riemann-Hilbert problem, while for other equations such as the Kadomtsev- 
Petviashvili II equation, this method is based on a certain generalization of the Riemann- 
Hilbert problem called the d (DEAR) problem. 



It is interesting that although equation (|T]) is an equation in two spatial variables, the 
Cauchy problem can be solved by a local Riemann-Hilbert problem. This is a consequence 
of the fact that the equation (|^) is a first order ODE in the variable x — ky. 

For integrable equations, there exist several different methods for constructing exact 
solutions. Such exact solutions for the Ward model in laboratory coordinates have been 
constructed in fl], . In particular. Ward constructed soliton solutions using the so-called 
dressing method 0. These solutions are obtained by assuming that M{x,y,t, k) has 
simple poles. In this case the corresponding solitons interact trivially, that is they pass 
through each other without any phase-shift. Recently, new soliton [0, ^ and soliton- 
antisoliton solutions ||^ were derived, by assuming that M{x, y, t, k) has double or higher 
order poles. The corresponding lumps interact nontrivially, namely they exhibit tt/N 
scattering between N initial solitons. 

The formalism presented in section 3 can also be used to obtain exact soliton solu- 



tions. In particular, it is shown in section 3 that if the assumption (|1^) is violated then 
M{x, y, t, k) still satisfies the Riemann-Hilbert problem (|1^) but now it is generally a 
meromorphic as opposed to a holomorphic function of k. The solitonic part of the solu- 
tion Q{x,y,t) is generated by the poles of M. The main advantage of this approach is 
that it can be used to establish the generic role played by the soliton solutions. Namely, it 
is well known |^ that the long time behaviour of the solution of a local Riemann-Hilbert 
problem of the type ([18|) where M is a meromorphic function of k, is dominated by the 
associated poles. Thus the long time behaviour of Q{x, y, t) with arbitrary decaying initial 
data Qo{x,y) is given by the multisoliton solution. 



2 THE CAUCHY PROBLEM WITHOUT SOLITONS 

In this section we prove Theorem 1. 

We first consider the direct problem, ie, we show that the spectral data S{x + ky, k) 
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are well defined in terms of the initial data Qo{x,y). Replacing Q{x,y,t) by Qo{x,y) in 
equation (|^) we find 

'-i^^-k?i^^-Q,^(.,y),(.,y,k) = 0. (23) 

Let fl{p, y, k) denote the x- Fourier transform of fi{x, y, k). Then equation (|23|) gives 



Jly-ipkjl- dxe Qox{x,y) ij{x,y,k) = 0. (24) 



R 



Equations ([l5|) and (0) are integrable forms of equation (^) with different initial values 



Under the small norm assumption (14), equations ([T5| ) and (16) are uniquely solvable in 



the space of bounded continuous functions f{x, y) such that / — / has a finite Li norm. 
Since the dependence on k of the kernel of the integral equations (|15|) and (p!6D is analytic, 
the functions fi'^{x,y, k) are analytic in k for ±Imk > 0. 

Equations (|1^) and (|16D can also be written in the form 

H^{x,y,k) = I+ dx'dy'G^{x-x',y-y',k)Qo^>{x',y') fi^{x',y',k), (25) 

where 

GHx,y,k) = —J^Jpdl^—^I, keC\ (26) 
We note that can be evaluated in closed form, 

G±(x, y, k) = ±^ + M ^e{x) - e{-x)) , (27) 

where S{y) and 9{x) denote the Dirac and the Heaviside functions, respectively. Indeed, 
writing 1/k = {ku — iki)/\k'^\ and using 

/ dx = 2m (sgnx) Oi-xh) e'^^^+'''\ a,b,p E R, (28) 

JR X + a + lb 

we find 

G^(x,y,k) = / dle'^yeixlki). (29) 

2-nk JR 

Recall that G^ corresponds to ki > 0; then in this case (sgnx) 6{xlki) = 6{x)6{l) — 
9{—x)9{—l), and the above equation becomes 

1 / /-oo , > rO 



i ( e{x) e{-x)\ 



2nk \y + iO y — iO 



(30) 
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Using 

^ = -T7rz%), (31) 
we find the expression for given by (^). Similarly for G~ . 



Using equation (|27|) it is straightforward to compute the large k behaviour of /x^: 

zttA; Ji?2 y — y 2k \j~oo Jx ) k^ 



for /c ^ 00. Thus 



/i = / ^ + '^P''' ^^oo- (33) 



Taking the complex conjugate of equation (plSf), letting p — > —p and using the fact 
that 

Qoii — ~Qoii = Qo22 5 Q012 ~ ~Qo2iy (34) 

we find 



(J'tiix, y, k) = /i22(a;, y, k), fiti{x, y, k) = -fii2{x, y, k), 



I^i2{x,y:k) = -fi2iix,y,k), /ij2(a;,?/, fc) = /iii(x,?/, /c). (35) 



Letting ^ = x + ky, rj = x — ky, equation (^) becomes 

2fc^-gox/i = 0. (36) 

Thus any two solutions of this equation are related by a matrix which is a function of 
X + ky and of k. Hence 

/i~(z, y, k) = fi'^{x, y, k) (/ — S{x + ky, k)) , k E R. (37) 



Equation (p7\} and the symmetry relations (pSf ), imply that / — S* is a Hermitian matrix. 
In particular, the determinant of / — S" is real. The determinant of equation ( p7|) yields 

det /X" = det /x+ det {I - S). (38) 



Taking the complex conjugate of this equation and using the sjnnmetry relations (pSf) , we 
find 

det /i+ = det /i" det (/ - 5) . (39) 
Equations ( ^8l) and (|39|) imply det(/ — 5) = ±1. However, equation (|33|) implies that 

det/x^ = l + 0(l/fc), k^oo. (40) 



Thus equation ( p8D implies det(/ — 5) = 1 + 0(1/ k) as k ^ oo, and since det (J — S') = ±1 
it follows that 

det(J -S) = l. (41) 

Equations (^S]) and ( ^TD imply 

det/i+ = det/i". (42) 



Since fi^ are analytic in C^, equations (^D|) and (^) define a local Riemann-Hilbert 
problem [|1^]. Its unique solution is 

det/i+ = det/i" = 1. (43) 

Evaluating equation (|37|) as y ^ — oo (keeping x + ky fixed) we find 

/ - ^ = f lim (/i+)-^) f lim (/i-)) , (44) 

\y->-oo } \y^-oo J 

which is equation (|l^). 

We now consider the inverse problem, ie, we show how to construct the solution of 
the Cauchy problem ([|), (|), starting from S{x + ky, k). Given S{x + ky + k'^t, k), we 
define M{x, y, t, k) as the solution of the Riemann-Hilbert problem (IT^) . In general, if the 



L2 norm with respect to of S and of || are sufficiently small, then the problem has a 
unique solution. However, in our particular case the solution exists without a small norm 
assumption. This is a consequence of the fact that / — 5* is a Hermitian matrix. Using 
this fact it can be shown (see for example [^) that the homogeneous problem, ie, the 
problem 

$- = $+(/ -5), keR, (45) 

$ = 0(^), A; ^00, (46) 

k 

has only the zero solution. 

Given M, we define Q{x,y,t) by equation (121]). A direct computation shows that if 
M"^ solves the Riemann-Hilbert problem (|TH), ie, if satisfies 



M+{x,y,t,k) = I + — dk — , k e R, 47 

2m JR k—{k + lO) 

and if Q{x, y, t) is defined by equation (|2lD then M+ satisfies equations (||) and (§). Hence 



Q satisfies equation (|l[) . Furthermore the investigation of the Riemann-Hilbert problem 



(p!8|) at if: = 0, implies that Q{x,y,0) = Qo{x,y). Also since J — S" is Hermitian, M+ and 
M~ have the proper symmetry properties (see equations (|35D), which in turn imply that 
Q{x,y,t) is a traceless ant i- Hermitian matrix. 

Remark It is important to note that equation (Rl) involves Qx and not Q. Because 



of this fact, equation ( |32D reduces to equation (|33|). This is to be contrasted with the 
Kadomtsev-Petviashvili equation, whose Lax pair involves Q. In that case equation (|32| ) 
simplifies to equation ( |55| ) only if JfiQ{x,y) dx = 0. This is the reason why it is usually 
assumed that the initial data of the Kadomtsev-Petviashvili equation satisfy the above 
condition. Without this assumption, the inverse spectral method is more complicated 



3 SOLITON SOLUTIONS 

In this section we show how the formalism of section 2 can be modified to include the 
soliton solutions. 

Since the matrix / — 5 is Hermitian of determinant one, it can be represented as 

where a is an arbitrary function of {x + ky + k'^t, k). 

Then equation (|T8]) becomes 

(Mf M^-) = (M+ M+) ^ 1 ) ' (49) 

where and M2 are 2-dimensional column vectors which are functions of {x,y,t,k). 
In particular 

Mf = M+ + aM+. (50) 
Equations ([T5|) and (16) are Fredholm integral equations of the second type; thus 



they may have homogeneous solutions. These homogeneous solutions which correspond 
to discrete eigenvalues are rather important since they give rise to solitons. We assume 
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that there exists a finite number of discrete eigenvalues and that they are all simple. Then 
Fredholm theory implies that admits the representation 

M+ = < + L k-k ' ^ ^ 

where mf{x,y,t,k) is analytic for k G and the vectors (f)i{x,y,t), I < I < N are 
homogenous solutions of the first column vector of equation (|I^). Following the arguments 
of it can be shown that 

(j)i{x,y,t) = ~ci {x + hy + /eft) M+ {x,y,t,ki), (52) 

where q is a scalar function of the argument indicated. Hence equation becomes 

. .J. +. .J. ^ci{x + ky + kft) M^{x,y,t,ki) 
M{{x,y,t,k) = mj{x,y,t,k) - }^ — . (53) 

1=1 '^i 

Substituting this equation into equation (^Dp solving the resulting Riemann-Hilbert prob- 
lem we find 

, ,^ Q(x + % + fcft)M+(x,i/,t,fc;) /l\ 1 f a{x + k'y + k'H)M+{x,y,t,k') 

(54) 

Let 

M+(x.y.a) = (i^f;/^';;j>). (55) 

In what follows, for simplicity of notion we suppress the x, y, t dependence. Using the 
notation (]55|) together with the symmetry relation (p5|), equation ([5^) becomes 



B{k) \ _( l\ f / A{h) \_J_f dk'ajk') ( A{k') \ 



-A{k) I \0 J t{k-ki\ B{ki) J 2m Jr k' - {k - iO) V B{k') 

Equation ( |5BD express the solution of the Riemann-Hilbert problem (^) in the case that 
solitons are included. 

PURE SOLITONS 

Soliton solutions correspond to a = 0. In this case evaluating equation (|56|) at = kj 
we find 

1=1 '^j '^i 
^ Q 



Aik,) = Y.^^B{h)- (57) 



1=1 kj h 
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The complex conjugate of these equations yields 



N 



B{k,) 



1-E 

1=1 ~ 



N 



Q 



1=1 "'J 



ki ki 



Equations ( 0) and ( pB] ) determine A{ki) and B{ki), I 



AT. 



(5^ 



Equation (|54D yields 



Thus using equation (^) we find 

N 



N 

E 



J f^^k-ki\ B{ki 



N 



Qii = -Q22 = E ^(^' 



1=1 



Q21 = -Ql2 = J2^lB{ki) 

1=1 



(59) 



(60) 



In summary, the A^-soliton solution is given by equations (|60|), where q = Ci{x + kiy + 
fcft) and A{ki), B{ki) are the solutions of the equations (|57D and (|58D. In the case of 
1-soliton equations (0) and (|58D yield 

1 Ci . 1 



A{k, 



1 ki-ki 



Bik, 



\ci\ 



{ki-kif 



Thus 



Qii = -Q 



1 



22 



1 _ fci - fci ' 



Q21 = -Q 



Cl 



12 



1 - 



\ci\ 



(61) 



(62) 



(fcl-fci)2 " (fcl-fcl)2 

Figure 1 represent a snapshot of the solution of equation (0) by taking ci = x + kiy + kft 
for /ci = i at time t = —3. 



4 THE WARD MODEL IN LABORATORY COORDINATES 

Here we study the Cauchy problem (H), (^) using the Lax pair (H), (|^). In the case that Qi 
Q2 are sufficient small, the inverse spectral method yields a solution of the Ward model 
in laboratory coordinates through the following construction. 

Theorem 2 Let Qi{x,y), Q2{x,y), x,y & R be 2 x 2 anti-Hermitian traceless matrices 
which are Schwartz functions and which satisfy the small norm conditions 

f QM,y) d^dy<^l, J = 1,2, (63) 

11 



Figure 1: 1-soliton solution of (|l|) at t = 



-3. 



where Qj is the Fourier transformation of Qj in the x variable. 

(i) Given Qi{x,y), Q2{x,y), define fi~^{x,y, k) , k G = {k E C : Irak > 0} and 
fi~{x, y,k), k ^ = {k E C : Irak < 0} as the 2x2 matrix valued functions which are 
the unique solutions of the linear integral equations 

47r\J0 J-oo J~oo Jy JJR \ k J 

(64) 

and 

47r\J-oo J-oo JO Jy JJR \ k J 

(65) 

(a) Given fi^ define the 2x2 matrix S{x + ^^^y, k), x,y,k E R, by 

(iv) Given S{x + ^-^y-, k) define the sectionally holomorphic function M{x,y,t, k) = 
M~^{x, y, t, k) for k G , M{x, y, t, k) = M~{x, y, t, k) for k G as the unique solution 
of the following 2x2 Riemann-Hilbert problem 

/ - 1 k"^ + 1 \ 

M- (x, y, t, k) = M+(x, y, t, k) 1 1 - S{x + + ^) 1 ' ^ ^ ^' (67) 
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detM=l, (68) 

M = I + 0{l/k), k^oo. (69) 

(v) Given M{x,y,t, k) define Q as 

Q,{x, y,t) = ^ f dk M+(x, y, t, k)S{x, y, t, k). (70) 



2m Jr 

Then Q solves equation ^ and Q{x,y,0) = Qi{x,y) , Qt{x,y,0) = Q2{x,y) . 

The proof of Theorem 2 is similar to section 2. 
Equations (|^) and (|65D can also be written in the form 

^i^ix,y,k)=I + J^^ dx'dy'G{x -x',y- y\ k) (q^,, + ~^''' ) ^^H^\y',k), (71) 



where 



or 



G±(i.,j/,*;) = ±^ + ^(»(x)-e(-x)), (73) 

Substituting this equation into equation ([7l|), it is straightforward to compute the large 
k behaviour of /x^, 



2T[k Jb2 y — y' 2A; \7-oo Jx J \ k J k"^ 

(74) 



for — > oo. Thus 

Qi {x,y) , 1 



/i = /-^^^ + 0(^), k^oo. (75) 



The corresponding soliton solutions of equation (||) can be derived following the 
method of section 3. 
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